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ABSTRACT 

We present a photometric study of the globular clusters (GCs) in the Virgo giant elliptical galaxy 
M86 based on Washington CT\ images. The colors of the GCs in M86 show a bimodal distribution 
with a blue peak at (C — T\) = 1.30 and a red peak at (C — T\) = 1.72. The spatial distribution of the 
red GCs is elongated similarly to that of the stellar halo, while that of the blue GCs is roughly circular. 
The radial number density profile of the blue GCs is more extended than that of the red GCs. The 
radial number density profile of the red GCs is consistent with the surface brightness profile of the M86 
stellar halo. The GC system has a negative radial color gradient, which is mainly due to the number 
ratio of the blue GCs to the red GCs increasing as galactocentric radius increase. The bright blue GCs 
in the outer region of M86 show a blue tilt: the brighter they are, the redder their mean colors get. 
These results are discussed in comparison with other Virgo giant elliptical galaxies. 
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1. INTRODUCTION 

Thousands of globular clusters (GCs) are found in 
a giant elliptical galaxy (gE) and they are found to be 
located from the center to the outer halo of their host 
galaxy (Lee 2003; Brodie & Strader 2006). Therefore, 
the GCs in gEs are a powerful tool to study the struc- 
ture and evolution of the GC system itself as well as 
their host galaxy. 

M86 (NGC 4406, VCC 881) is a famous gE located 
close to the center region of the Virgo galaxy cluster 
in the sky and is known to be infalling to the Virgo 
center with the relative velocity of about 1300 km s _1 . 
It is also known early to be abundant with GCs (Hanes 
1977). There are several photometric studies for the 
GCs in M86 in the literature. Cohen (1988) presented 
gri photometry of GCs at R < 7' in M86 obtained at 
the Hale 5 m telescope. She found that the GC system 
of M86 is similarly extended as the stellar light, and 
that there is no detectable radial color gradient. 

From HST/WFPC2 VI observation of the central 
region of M86, Neilsen & Tsvetanov (1999) reported 
that the color distribution of the GCs in M86 shows 
a single peak. However, Kundu & Whitmore (2001) 
and Larsen et al. (2001) showed later, using the same 
data as used in HST/WFPC2 by Neilsen & Tsvetanov 
(1999), that the color distribution of the GCs in M86 
is bimodal. Later Peng et al. (2006) confirmed, from 
HST/ACS gz observation, that the color distribution 
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of M86 GCs is bimodal. 

Jordan et al. (2007) and Villegas et al. (2010) de- 
rived a luminosity function of M86 GCs from the 
same HST/ACS images, and obtained a similar Gaus- 
sian peak value at g-band, /j, g — 23.950 ± 0.097 and 
23.887 ± 0.087, respectively. Peng et al. (2008) es- 
timated the specific frequency of M86 GCs, Sn = 
2.57 ± 0.12, from the HST/ACS images. On the other 
hand, Rhode & Zepf (2004) investigated M86 GCs us- 
ing the wide field (36' x 36') BVR images obtained 
using the KPNO 4 m Mayall telescope. They found 
again that M86 GCs show a bimodal color distribution 
and that they show a modest negative color gradient 
with galactocentric radius. They derived a specific fre- 
quency of Sn = 3. 5 ±0.5, which is larger than the value 
derived from the central region by Peng et al. (2008). 

However, previous studies based on HST data cov- 
ered only a small central region of M86 and those based 
on ground-based data did not investigate the prop- 
erties of sub-populations in the GC system and the 
spatial distribution of the GCs in M86. Here we in- 
vestigate the detailed properties of the M86 GC sys- 
tem using Washington photometry derived from deep 
and wide CCD images. The Washington filter sys- 
tem with a wide bandwidth is known to be very sen- 
sitive to measuring the metallicity of the GCs (Geisler 
& Forte 1990) so that it is useful for studying the 
properties of sub-populations of the GCs. We adopt 
a distance to M86, derived from the surface bright- 
ness fluctuation method by Mei et al. (2007), 16.9 Mpc 
((m — Al)o = 31.13 ± 0.07). At this distance one arc- 
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sec corresponds to a linear scale of 81 pc. The basic 
information of M86 is listed in Table 1. 

This paper is composed as follows. Section 2 de- 
scribes observation and data reduction. In §3, we 
present the color-magnitude diagram and color distri- 
bution of the M86 GCs. We investigate the spatial 
distribution of the GC system as well as the radial vari- 
ation of number densities and colors of the GCs. In §4, 
we discuss our results and their implication in compar- 
ison with other Virgo gE studies. Primary results are 
summarized in §5. 

2. OBSERVATION AND DATA REDUCTION 

2.1 Observations 

The images of M86 were obtained with the 2048 x 
2048 pixels CCD camera (the pixel scale is 0.47 arcsec 
pixel -1 ) at the KPNO 4 m telescope on the photomet- 
ric night of April 9, 1997. These images were taken 
with the Washington C and Kron-Cousins R filters. 
The R filter was used instead of the Washington T\ fil- 
ter, as done in most recent Washington studies (Geisler 
1996; Lee et al. 1998; Harris et al. 2004; Dirsch et al. 
2005; Bassino et al. 2006; Forte et al. 2007; Lee et al. 
2008), because the sensitivity of the Kron-Cousins R 
filter is about three times larger than that of the T\ 
filter, while the effective wavelengths of two filters are 
very similar as R = T\ + 0.003 - 0.017(C - Ti) with 
rms=0.02 (Geisler 1996). Hereafter we call this R filter 
as T\ filter in the rest of this paper. 

Table 2 shows the observation log. The size of the 
field of view is 16' x 16', and the seeing range is from 
1".0 to 1".2. Exposure times are 60 s and 4 x 1500 s 
for C, and 30 s and 2 x 1000 s for T x . Several Washing- 
ton standard fields in Geisler (1996) were also observed 
during the observing run. 

2.2 Point Source Photometry 

Using the IRAF software, each frame was bias- 
subtracted and flat-fielded. Then the long exposure im- 
ages of two Ti-bands and four C-bands were combined 
with average and median procedure, respectively. Fig- 
ure 1 displays a T\ image of M86 taken with the short 
exposure. 

The GCs at the M86 distance appear as point 
sources in the KPNO images so that we derived pho- 
tometry of the point sources from the images as fol- 
lows. First, we derived isophotal model images of M86 
with ELLIPSE task in IRAF/STSDAS. Then we sub- 
tracted these images from the original images to detect 
better point sources in the images. Next we carried 
out the point-spread function (PSF) photometry on the 
subtracted images using the DAOPHOT/ALLFRAME 
(Stetson 1994). Then we derived the aperture cor- 
rection using the photometry of several isolated bright 
stars and applied it to the PSF-fitting magnitudes: 
0.013 ±0.006 and 0.019 ±0.004 for short and long expo- 
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Fig. 1. — A grey scale map of Ti image of the M86 field 
taken with KPNO 4 m telescope. The field of view is 16' x 
16'. The boxes represent the HST/ACS fields. The galaxy 
in the north is NGC 4402, which is an edge-on galaxy of Sb 
type. 



sure C images and -0.025 ±0.008 and 0.022 ±0.003 for 
short and long exposure T\ images. Finally, we trans- 
formed the instrumental magnitudes of the sources onto 
the standard system using the transformation equa- 
tions (April 9) given in Section 2.3. by Lee et al. 
(2008). We selected the point sources among the de- 
tected sources using the morphological classifier r_2 
moment (Kron 1980). We considered the sources with 
r_ 2 < 1.27 as the point sources. Figure 2 displays the 
mean errors of T\ and (C — Ti) versus T\ magnitude, 
and Table 3 lists the CT\ photometry catalog of only 
4351 point sources with error (C — Ti) < 0.3 measured 
in the KPNO images. 

The central region at R < 1 arcmin in the long ex- 
posure images was saturated and the short exposure 
images are too shallow. Therefore, we used HST pho- 
tometry available in Jordan et al. (2009) for the anal- 
ysis of the central region of M86. Jordan et al. (2009) 
provides a photometric catalog of GCs in 100 early- 
type galaxies in Virgo including M86, derived from 
the HST/ACS Virgo Cluster Survey (ACSVCS). Their 
catalog includes g and z magnitudes, and the proba- 
bility for GC classification (pgc)- GCs appear as ex- 
tended objects in the ACS images and the sources with 
Pgc > 0.5 are considered as GC candidates. 

We derived the transformation between (C — Ti) 
colors and (<? — z) colors for the bright objects (with 
Ti < 22.0 mag, error (C - Ti) < 0.08 and 1.1 < 
(C — Ti) < 2.1), using the common objects between 
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Table 1. 

Basic information of M86 



Parameters 




Values 


References 


R.A., Decl. (J2000) 


12 


h 26 m 11.743 s , +12° 56' 46.40" 


1 


Total magnitudes 


V T 


= 8.90 ± 0.05, B T = 9.83 ± 0.05 


2 


Foreground reddening, E(B — V) 




0.030 


3 


Distance, d 


16.8f 


i Mpc ((to - M) = 31.13 ± 0.07) 


4 


Systemic radial velocity, v p 




-244 ± 5 km s" 1 


5 


Effective radius, R e ff 


3.59 arcmin (C), 3.14 arcmin (Ti) 


6 


Effective ellipticity, e e // 




0.36 (C), 0.33 (Ti) 


6 


Effective P. A., 6 e// 




119 deg (C,Ti) 


6 


Standard radius, R25 


4.8 


'.8 arcmin (C), 9.64 arcmin (Ti) 


6 


Standard ellipticity, £25 




0.40 (C,Ti) 


6 


Standard P. A., 25 




124 deg (C,Ti) 


6 



References. — (1) NASA Extragalactic Database; (2) de Vaucouleurs et al. (1991); (3) 
Schlegel et al. (1998); (4) Mci et al. (2007); (5) Smith et al. (2000); (6) This study. 



Table 2. 

Observation log 



Target 


Filter 


T(cxp) 


Airmass 


Seeing 


Date 






(sec) 




(") 


(UT) 


M86 


C 


60 


1.12 


1.13 


1997 Apr 9 


M86 


R 


30 


1.06 


0.99 


1997 Apr 9 


M86 


C 


4 x 1500 


1.06 


1.13 


1997 Apr 9 


M86 


R 


2 x 1000 


1.21 


1.22 


1997 Apr 9 
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Fig. 2. — Mean photometric errors of Ti magnitude and 
(C — Ti) color for the point sources obtained from the long 
exposure images. The circles with error bars, triangles, and 
squares represent the mean errors of all point sources, the 
blue GCs (1.0 < (C-Ti) < 1.55), and the red GCs (1.55 < 
(C-Ti) < 2.1), respectively. 



the KPNO results and the ACSVCS catalog. Figure 3 
shows relations for colors and magnitudes. There are 
several outliers, which are due to being located near 
saturated sky levels or having nearby sources. We fitted 
the data after excluding these outliers. Linear fitting 
to the data yields [g - z) = 0.773(C - T x ) - 0.118 with 
rms=0.068, and g = 0.998Ti + 0.911 with rms=0.109. 

2.3 Completeness of the Photometry 

We estimated the completeness of the KPNO pho- 
tometry using DAOPHOT/ADDSTAR, which is de- 
signed for the artificial star experiment. First we 
generated a set of artificial stars of which the color- 
magnitude diagrams and luminosity functions are sim- 
ilar to those of the detected point sources, using the 
PSFs derived from the long exposure images. Next 
we added them to the real image avoiding the posi- 
tion of detected real sources. We obtained 50 artificial 
images each of which includes 1200 artificial stars for 
each filter so that the total number of added artificial 
stars is 60,000. Then we applied the same photome- 
try procedure to the artificial images as used for the 
real image. From these results, we estimated the com- 
pleteness, which is the number ratio of the recovered 
artificial stars to the added artificial stars. 

Figure 4 displays the completeness of C and T\ mag- 
nitude measured for the point sources. This figure 
shows the completeness for the entire region and four 
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Fig. 3. — Comparison of (C — Ti) with (g — z) (a) and Ti 
magnitude with g magnitude (b). The solid lines represent 
the linear least-squares fits for error (C — Ti) < 0.08 and 
Ti < 22.0, except for the outliers indicated by the crosses. 



radial bins: 1' < R < 3', 3' < R < 5', 5' < R < 7', 
and 7' < R < 9'. The 90 % completeness is located at 
Ti = 23.8 (C = 25.1) for the outer region at R > 3' 
and T\ = 23.5 (C = 24.9) for the inner region at 
1' < R < 3'. The photometric limit levels with 50 
% completeness are Ti = 24.7 (C = 25.9) for the outer 
region at R > 3' and T x = 24.4 (C = 25.6) for the 
inner region at 1' < R < 3'. The completeness varies 
little, depending on galactocentric radius for the outer 
region at R > 3', while it is somewhat lower for the 
inner region at R < 3' compared with the outer region. 

3. RESULTS 

3.1 Color-Magnitude Diagram 

Figure 5 displays the color-magnitude diagrams 
(CMDs) of the point sources measured from the KPNO 
images of M86 as well as that from the HST/ACS pho- 
tometry (for R < 1'.5). We plotted the CMDs for the 
M86 region (1'.5 < R < 10') and the background region 
(10' < R < 12'). For the central region at R < 1'.5 we 
plotted only the HST photometry. We transformed the 
(g — z) colors and g magnitudes of the objects given in 
the ACSVCS catalog into (C— Ti) colors and T\ magni- 
tudes, respectively, using the transformation equations 
derived in the previous section. 
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Table 3. 

A catalog of Washington photometry for the point sources in M86 a 



ID 


R.A. 


Decl. 


Ti±error 


(C - Ti)±error 




(J2000) 


(J2000) 


(mag) 


(mag) 


1 


12:26:21.75 


12:52:40.7 


17.746 ± 0.011 


2.883 ±0.033 


2 


12:26:03.70 


12:59:57.5 


17.888 ±0.023 


2.750 ±0.030 


4 


12:26:13.36 


13:05:11.7 


17.992 ±0.010 


0.777 ±0.015 


5 


12:26:28.38 


13:00:52.8 


18.045 ±0.014 


0.792 ±0.025 


6 


12:26:17.03 


13:06:00.4 


18.195 ± 0.011 


2.163 ±0.022 



a The complete version of this table is in the electronic edition of the Journal. 
The printed edition contains only a sample. 




22 24 26 28 

C 




T, 

Fig. 4. — Completeness of C magnitude (a) and T\ mag- 
nitude (b). The solid lines represent the results for the 
entire region. The short dashed, long dashed, dot-short 
dashed, and dot-long dashed lines indicate the complete- 
ness at 1' < R < 3', 3' < R < 5', 5' < R < 7', and 
7' < R < 9', respectively. 



Three distinct features are seen in Figure 5: (a) GCs 
of M86 in the broad vertical feature in the color range 
of 1.0 < (C - Ti) < 2.1, (b) a small number of bright 
foreground stars at (C - Ti) < 1.0 and (C - Ti) > 2.5, 
and (c) faint blue unresolved background galaxies with 
T\ > 23 mag. We selected the bright point sources 
with 1.0 < (C - Ti) < 2.1 and 19 < T x < 23.08 for 
the analysis of the M86 GCs, considering the following 
points: (1) the color range of the known GCs in other 
galaxies is (C - T{) « 0.9 to 2.1 (Geisler 1996; Lee 
et al. 1998; Dirsch et al. 2005; Forte et al. 2007; Lee 
et al. 2008); (2) the peak luminosity of the M86 GCs is 
T\ = 23.08, which was converted from the peak g mag- 
nitude (/i g = 23.95±0.10) given by Jordan et al. (2007); 
(3) the mean photometric error and the completeness 
for Ti < 23.08 are smaller than error (Ti) = 0.05 and 
larger than 90%, respectively; and (4) the contamina- 
tion due to background galaxies is estimated to be very 
small for Ti < 23.08. 

3.2 Color Distribution of the Globular Clus- 
ters 

Figure 6 shows the color distribution of the point 
sources with 19 < Ti < 23.08 at 1'.5 < R < 10' 
in the KPNO images. The color distribution of the 
background objects with the same magnitude range at 
10' < R < 12' is also shown. 

We carried out the KMM test, which estimates the 
statistical significance of bimodality in astronomical 
datasets (Ashman et al. 1994), to investigate the bi- 
modality of the (C — Ti) color distribution of the M86 
GCs. KMM test of the data shows that the probabil- 
ity for bimodal color distribution is higher than 99.9 
%. The color distribution is also well described as the 
best-fit double Gaussian curves: a blue component with 
center at (C — T\) = 1.30 and width a = 0.15 and a red 
component with center at (C — Ti) = 1.72 and width 
a = 0.15. The minimum between the two components 
is found to be at (C — Ti) w 1.55. Thus we divided the 
entire sample of the M86 GCs into two sub-populations: 
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Fig. 5. — Color-magnitude diagrams of the sources in M86 given in ACSVCS catalog (a) and the point sources in 
M86 measured from the KPNO images (b and c). The boxes indicate the color (1.0 < (C — 21) < 2.1) and magnitude 
(19 < T\ < 23.08) criteria for selecting the M86 GCs. The dotted and crosses in panel (a) represent the values converted 
from g and (g — z) at R < 1'.5 and l'.5 < R < 2'. 4, respectively. 
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Fig. 6. — Color distribution of the M86 GCs with 
19 < Ti < 23.08 and l'.5 < R < 10'. The solid line and 
dashed lines represent the double Gaussian fit derived from 
the KMM test. The hashed histogram is the color distri- 
bution of the point sources within the background region 
(10' < R < 12'). The vertical dotted lines indicate the 
color boundaries for the blue GCs (1.0 < (C - Ti) < 1.55) 
and the red GCs (1.55 < (C - Ti) < 2.1). 



the blue GCs (BGC) with 1.0 < (C - Ti) < 1.55 and 
the red GCs (RGC) with 1.55 < (C - Ti) < 2.1. This 
boundary color is consistent with the previous study of 
M86 (Rhode & Zepf 2004), but is bluer than that used 



for M60 (Lee et al. 2008) and M49 (Lee et al. 1998) 
and slightly redder than that used for M87 (Cote et al. 
2001). Figure 7 displays the radial variation of the 
color distribution of the M86 GCs. The number ra- 
tio of the blue GCs and the red GCs increases as the 
galactocentric radius increases. 

3.3 Surface Photometry of M86 

We derived surface photometry of M86 from the 
KPNO images using the IRAF/ELLIPSE task. Af- 
ter masking out the bright foreground stars and back- 
ground galaxies and using the mode values of the re- 
gions at R J> 10' in the north-east corner as the back- 
ground levels, we carried out ellipse fitting to the C 
and Ti images. The errors for the surface brightness 
magnitudes are given by ELLIPSE. The long exposure 
images are saturated at the central region of M86, while 
the short exposure images are not. Thus we used the 
values from the short exposure images at R < 1', the 
average values from the short and long exposure images 
at 1' < R < 3', and the values from the long exposure 
images at R > 3'. 

Figure 8 shows the radial profiles of the surface 
brightness magnitudes, (C — Ti) color, ellipticity, and 
position angle (P.A.) of M86 as a function of major ra- 
dius R ma j ■ The surface brightness profiles of the C and 
T\ magnitudes are consistent with that of the B mag- 
nitude given by Caon et al. (1990). The surface bright- 
ness profiles are roughly described by a de Vaucouleurs 
law. These surface brightness profiles are fitted with 
a de Vaucouleurs i? 1 / 4 law for the range of 2" < 
Rmaj < 400", using linear least-squares fitting: fi(C) = 
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Fig. 7. — Radial color distribution of the GCs in M86. 
The histograms in each panel and the dot-dashed line with 
the circles in panel (a) represent the color distributions of 
the GCs obtained from KPNO and from ACSVCS catalog, 
respectively. The vertical dotted lines indicate the color 
boundaries for the sub-populations. 



2.173(±0.032)i? 1 / 4 + 16.009(±0.086) with rms=0.219, 
and n(Ti) = 2.247(±0.033)i? 1 / 4 + 13.981(±0.089) with 
rms=0.226. The effective radius (R e ff) and standard 
radius (R25) for C-band are derived to be 3'. 59 and 
4'. 88 (3'. 14 and 9'. 64 for Ti), corresponding to linear 
sizes of 17.54 kpc and 23.85 kpc, respectively. 

This value for the effective radius is much smaller 
than that derived from the small field HST/ACS im- 
ages by Ferrarese et al. (2006) (R eS = 6'.86 ~ 33.52 
kpc at <?-band) but similar to that derived from wide 
field images by Kormendy et al. (2009) (R c s = 3'. 38 ~ 
16.52 kpc at F-band). The surface brightness pro- 
files of C and Ti magnitudes in this study as well as 
the B magnitude given by Caon et al. (1990) show a 
break at R m aj ~ 135", resulting in flatter profiles in 
the outer region. Thus we carried out double linear 
fits with i? 1 / 4 law at 10" < R ma] < 135": (i(C) = 
2.240(±0.016)i? 1 / 4 + 16.027(±0.042) with rms=0.038 



and /i(Ti) = 2.304(±0.021)i? 1 / 4 + 14.030(±0.052) with 
rms=0.047, and 135" < R maj < 560": fi(C) = 
1.706(±0.042)i? 1 / 4 + 17.684(±0.170) with rms=0.063 
and /i(Ti) = 1.860(±0.056)i? 1 / 4 +15.338(±0.227) with 
rms=0.084. 

The (C — Ti) color of the M86 nucleus is very red 
with (C — Ti) R3 2.1. It becomes rapidly bluer in the 
inner region (R m aj < 50") as the galactocentric radius 
increases, but changes slowly in the outer region. This 
trend is consistent with the radial variation of (U — R) 
colors given by Peletier et al. (1990) as shown in Fig- 
ure 8 (b), but is somewhat different from the (B — R) 
colors given by Peletier et al. (1990). This may be re- 
lated with the existence of dust close to the M86 cen- 
ter, as seen in the far-infrared images (Stickel et al. 
2003; Gomez et al. 2010). The color at 2" < R maj < 
300" is approximately fitted by the log-linear relation 
(C - Ti) = -0.105(±0.003)logi? + 1.999(±0.005) with 
rms=0.014. The ellipticity increases rapidly from 0.15 
in the central region to 0.25 at R ma j ~ 40", staying 
almost constant at 40" < R m aj < 135". Then it in- 
creases again to ~ 0.4 in the outer region. Thus the 
isophotes of the outer region are more elongated than 
the central region. The values of the ellipticity at the 
effective radius and standard radius are 0.35 and 0.40, 
respectively. The value for P.A. is constant with 120±6 
deg. 

3.4 Spatial Distribution of the Globular Clus- 
ters 

We investigated the spatial structure of the M86 GC 
system. Figure 9 displays the spatial distributions of 
the GCs with 19 < Ti < 23.08 mag (all GCs, blue 
GCs, and red GCs) as well as M86 stellar light. We 
also overlay number density contours of the GCs in 
this figure. T% image of M86 with isophotal contours is 
displayed in the panel (b) for comparison with the GC 
system. 

Several features are noted in this figure. First, all 
the GCs have a circular spatial distribution, and show 
a central concentration. Second, the blue GCs have a 
roughly circular spatial distribution, extended farther 
than the red GCs. Third, the red GCs have a spatial 
distribution somewhat elongated along the southeast- 
northwest direction, which is consistent with the dis- 
tribution of M86 stellar halo. Using the dispersion 
ellipse method (Trumpler & Weaver 1953; Hwang & 
Lee 2008; Strader et al. 2011), we measured e (ellip- 
ticity) = 0.04 ± 0.02, 0.02 ± 0.03, and 0.16 ± 0.05, and 
P.A.=122 ± 24, 42 ± 43, and 131 ± 10 deg for all the 
GCs, the blue GCs, and the red GCs at R < 5', re- 
spectively. The value for the ellipticity of the red GCs 
is eight times larger than that of the blue GCs. The 
P.A.'s of all the GCs and the red GCs are similar to 
that of M86 stellar light. 
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Fig. 8. — Surface photometry of M86. Rmaj is galactocentric distance along the major axis of M86 stellar light, (a) 
Surface brightness magnitudes for Ti and C . The solid lines represent the linear least-squares fits with R 1 ^ law for 
2" < R ma j < 400". The dot-dashed lines represent the double linear fits with R 1/4 law for 10" < R ma j < 135" and 
135" < Rmaj < 560". The crosses indicate the surface brightness profile for B magnitude given by Caon et al. (1990). (b) 
Surface color profiles. The solid line indicates the power law fit of (C — Ti) and log(ii) for 2" < Rmaj < 300". The circles 
and crosses represent the colors from this study and from Peletier et al. (1990), respectively, (c) Ellipticity versus Rmaj- 
(d) Position angle (P. A.) versus Rmaj- 



3.5 Radial Number Density Profiles of the 
Globular Clusters 

We derived the radial number density profiles of the 
GCs with 19 < Ti < 23.08 mag, using the KPNO re- 
sults for the outer region at R > 1'.5 and the ACSVCS 
catalog for the central region at R < 1'.5. Here we 
masked out the region with 5' x 1'.5 centered on NGC 
4402 to remove any contamination due to this galaxy. 
First, we applied the correction for the incomplete- 
ness of photometry using the results shown in Figure 4. 
Rhode & Zepf (2004) found that the radial number den- 
sity of the GCs in M86 becomes zero at R > 17'. The 
field of view of our KPNO images is not wide enough to 
cover the entire region of M86. Here we adopted, as a 
background number density, a half of the number den- 
sity of the point sources at 10' < R < 12': 0.915±0.279 
arcmin -2 for all the GCs, 0.767±0.255 arcmin -2 for the 
blue GCs, and 0.149 ±0.113 arcmin- 2 for the red GCs. 
We subtracted these values from the original counts. 
Then we derived radial number densities, considering 



the effect of the GC luminosity function (by multiply- 
ing by a factor of two) . The radial number densities of 
the GCs thus derived are listed in Table 4. 

Figure 10 displays the radial number density profiles 
of the GCs versus i? 1 / 4 in (a) and R in (b) and (c). 
We also plot the radial number density profile of M86 
GCs given by Rhode & Zepf (2004) and the C-band 
surface brightness profile of M86 in (a) for comparison. 
Several notable features are seen as follows. First, the 
density profile of all the GCs is more extended than 
the galaxy stellar light. The density profile of all the 
GCs derived in this study agrees well with that given by 
Rhode & Zepf (2004). Second, the density profile of the 
blue GCs decreases more slowly as the galactocentric 
distance increases than that of the red GCs. Third, 
the density profile of the red GCs is consistent with 
the surface brightness profile of the galaxy stellar light. 

Fourth, the density profiles of the GCs at 1' < 
R < 10' are approximately fitted either by a de Vau- 
couleurs R 1 ^ law: logn GC = -1.683(±0.104) J R 1 / 4 + 
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Fig. 9. — Spatial distributions and number density contours of M86 GCs. The plus signs indicate the center of the galaxy, 
(a) All GCs. The contour levels are 2.037, 5.093, 8.148, and 12.223 GCs arcmin -2 . (b) The grey map with isophotal contours 
in the T\ image of M86. (c) Blue GCs. The contour levels are 0.6 times of those for all the GCs. (d) Red GCs. The contour 
levels are 0.4 times of those for all the GCs. 



3.379(±0.154) for all the GCs, log n GC = -1.534(±0.082) 
i? 1 / 4 + 3.006(±0.120) for the blue GCs, and logn GG = 
-2.048(±0.182) J R 1 / 4 + 3.382(±0.268) for the red GCs; 
or by a power law: logn GG = — 1.343(±0.116) logi? + 
1.783(±0.080)forallthe GCs, log n GG = -1.226(±0.094) 
logi?+ 1.552(±0.064) for the blue GCs, and logn GG - 
-1.628(±0.185) logi?+ 1.434(±0.127) for the red GCs. 
In the case of all the GCs, the values for the slope agree 
well with those given by Rhode & Zepf (2004). Effec- 
tive radii of all the GCs, the blue GCs, and the red 
GCs are derived from the de Vaucouleurs i? 1 / 4 law fits: 
15'.33 (= 74.50 kpc), 22'. 20 (= 107.89 kpc), and 6'.99 
(= 33.97 kpc), respectively. Thus the effective radius 
of the red GCs is much smaller than that of the blue 
GCs, showing the central concentration of the red GCs 
is much stronger than that of the blue GCs. The den- 



sity profiles of the GCs show a break at R « 4', result- 
ing steeper profiles in the outer region. The outer parts 
at R > 4' are fitted well by a power law: logn GG — 
-1.915(±0.085) logi? + 2.255(±0.073) for all the GCs, 
logn GG = -1.672(±0. 101) logi? + 1.923(±0.086) for 
the blue GCs, and logn GG = -2.666(±0. 109) logi? + 
2.276(±0.093) for the red GCs. The inner parts at 
i? < 4' are fitted well by a power law: logn GG = 
-0.569(±0.094) logi?+ 1.496(±0.045) for all the GCs, 
logn GG = -0.596(±0. 100) logi? + 1.322(±0.047) for 
the blue GCs, and logn GG = -0.516(±0.089) logi? + 
1.014(±0.042) for the red GCs. In the case of the outer 
parts, the radial number density profile of the red GCs 
is much steeper than that of the blue GCs, while in 
the case of the inner parts, the sub-populations have 
similar slopes. 
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Fig. 10. — Radial number density profiles of M86 GCs. 
The filled and open points are from HST and KPNO, re- 
spectively. The circles, triangles, and squares represent all 
the GCs, the blue GCs, and the red GCs, respectively, (a) 
Number density versus R 1 ^ 4 . The small filled circles in- 
dicate the number density profile given by Rhode & Zepf 
(2004). The solid lines represent the fits for a de Vau- 
couleurs law for 1' < R < 10'. The dot-dashed line in- 
dicates the C surface brightness profile of the M86 stellar 
halo, (b) Number density versus R. The solid lines rep- 
resent the fits for a power law. The dotted lines indicate 
the two-component power law fits for O'.l < R < 4' and 
4' < R < 12'. (c) Number density versus R. The solid lines 
represent King model fits for R < 4', which are represented 
by the dotted lines for R > 4'. 



Finally, the density profiles of the GCs are some- 
what flat in the central region at R 2'. This flat ra- 
dial profile, which may indicate dynamical relaxation 
of the GC system (Cote et al. 1998) or the relic of the 
GC formation epoch (Harris et al. 1998; Kundu et al. 
1999), can be approximately fitted by a King model 
(King 1962). Figure 10 (c) displays the results of a 
King model fitting for R < 4'. We derived core radius 
r c = 1'.82 and concentration parameter c = 5.13 for all 
the GCs, r c — 1'.71 and c = 5.34 for the blue GCs, and 
r c = 2'. 07 and c = 5.68 for the red GCs. It is noted 



that in the case of the blue GCs the radial density pro- 
file of the outer region at R > 4' shows an excess over 
the King model. 

3.6 Radial Variation of Globular Cluster Col- 
ors 

We investigated the radial variation of the color of 
the M86 GCs. Figure 11 displays the (C-Ti) col- 
ors of the GCs with 19 < Ti < 23.08 mag as a 
function of galactocentric distance derived from the 
KPNO images. The mean color of all the GCs in 
each radial bin shows a radial gradient, while those 
of the blue GC and the red GC show little ra- 
dial gradients. The linear least-squares fitting at 
R < 10' yields (C - Ti) = -0.016(±0.003)-R + 
1.512(±0.018) with rms=0.026 for all the GCs, (C - 
Ti) = -0.008(±0.002)i2+1.341 (±0.012) with rms=0.017 
for the blue GCs, and (C - Ti) = 0.002(±0.005)i? + 
1.743(±0.028) with rms=0.040 for the red GCs. The 
color gradient of all the GCs is due to the variation of 
the number ratio of the blue GCs to the red GCs as 
shown in Figure 7. This negative color gradient for all 
the GCs is consistent with the previous studies based 
on (B-R) color (Rhode & Zepf 2004). The mean color 
of the M86 stellar light is much closer to that of the red 
GCs than to that of the blue GCs, as often seen in other 
gEs such as M49 and M60 (Lee et al. 1998, 2008). 

4. DISCUSSION 

We compared the photometric results of the M86 GC 
system with those of GC systems in other two gEs in 
Virgo (M60 and M49). M60, which is the third bright- 
est gE in Virgo and located about 3 deg from the Virgo 
center, and M49, which is the brightest gE in Virgo and 
located about 4 deg from the Virgo center. Previous 
studies for the GC systems of M60 (Lee et al. 2008) 
and M49 (Geisler et al. 1996; Lee et al. 1998) used the 
same instruments and similar photometric techniques 
as this study for the M86 GC system. For comparison 
we used the same criteria for the GC sample in three 
galaxies, GCs with R < T and 19 < Ti < 23. 

Figure 12 displays the color distributions of the GCs 
in M86, M60, and M49. The GCs in these galaxies 
show similarly bimodal color distributions, but there 
are slight differences among these galaxies. The peak 
colors for M86 ((C-Ti) = 1.31 and 1.74) are similar to 
those for M49 ((C-Ti) = 1-30 and 1.79), but are about 
0.1 mag bluer than those for M60 ((C-Ti) = 1-37 and 
1.87). Applying this color difference to the relation 
between (C— Ti ) color and metallicity (Lee et al. 2008), 
this implies that the M86 and the M49 GCs are on 
average about 0.2 dex more metal-poor than M60 GCs. 
These results are approximately consistent with those 
by Peng et al. (2006) who covered the central regions 
with R <J 1'.5. They reported that the peak colors for 
these three galaxies are {g — z) = 0.98 and 1.33 for M86, 
( g - z ) = 0.97 and 1.42 for M49, and (g - z) = 0.98 
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Fig. 11. — (C - Ti) color of M86 GCs as a function of 
galactocentric distance. The circles, triangles, and squares 
represent all the GCs, the blue GCs, and the red GCs, re- 
spectively. The filled and open points represent mean and 
median color of the GCs in each radial bin, respectively. 
The solid lines are the linear least-squares fits about the 
mean color of the sub-populations. The curved line indi- 
cates (C — Ti) surface color profile of M86 stellar halo. 



and 1.45 for M60. 

There are significant differences among the number 
ratios (N(BGC)/N(RGC)) of the blue GCs to the red 
GCs in these galaxies: 2.37 ± 0.26 for M86, 1.38 ± 0.11 
for M60, and 0.97 ± 0.07 for M49. Thus the fraction of 



the blue GCs for M86, the most faint among these three 
galaxies, is the largest. This is consistent with previous 
finding that the fraction of the blue GCs decreases as 
their host galaxy gets brighter (Peng et al. 2006). On 
the other hand, the number ratio of the blue GCs to 
the red GCs in each galaxy is a little bit larger than the 
value by Peng et al. (2006) ( 2.33±1.15 for M86, 0.75± 
0.10 for M60, and 0.69 ±0.09 for M49) and by Faifer et 
al. (2011) (1.07 ± 0.10 for M60). The radial coverages 
used by Peng et al. (2006) and Faifer et al. (2011) are 
R < 1'.5 and R <; 5', respectively, so that their fields 
are on average closer to the galaxy center than that 
used in this study (1' < R < 7'). These shows that 
this number ratio decreases as the galactocentric radius 
decreases (see Section 3.2 for detail). 

Recent studies for the GCs in gEs based on HST /ACS 
data (Harris et al. 2006; Strader et al. 2006; Mieske 
et al. 2006; Peng et al. 2009; Mieske et al. 2010) and 
wide-field images (Forte et al. 2007; Lee et al. 2008; 
Harris 2009a), found a so-called 'blue tilt' (or a color- 
magnitude relation): the brighter the bright blue GCs 
are, the redder their colors get. Several mechanisms 
were proposed to explain the origin of the blue tilt 
(Mieske et al. 2006; Bekki et al. 2007; Strader & Smith 
2008; Harris 2009b; Blakeslee et al. 2010). The self- 
enrichment process among these appears to be a main 
driver (Mieske et al. 2006; Strader & Smith 2008). We 
investigated any existence of this blue tilt in the M86 
GCs. Figure 13 displays the CMDs of M86, M60, and 
M49 derived from the CT\ images for the outer regions 
at 1'.5 < R < 7'. The blue and red peak values are 
determined from the KMM test with the colors in each 
magnitude bin. In Figure 13 the blue tilt is seen for the 
bright blue GCs in M86 as well as other two galaxies. 
However, the number of galaxies is only three so that 
it needs to study more galaxies to address the slope 
difference of the blue tilt. The existence of the blue tilt 
in M86 as well as M49 and M60 is consistent with the 
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Fig. 12. — Comparison of the color distribution between 
the GCs in M86 and the GCs in other Virgo gEs. (a) M86 
from this study, (b) M60 from Lee et al. (2008). (c) M49 
from Geisler (1996) and Lee et al. (1998). The solid and 
dashed lines represent the double Gaussian fits. 



result for high mass galaxies in the AC SVC S given by 
Mieske et al. (2010). 

The spatial distribution of the GCs in three galaxies 
share common features as follows: (a) the radial num- 
ber density profile of the red GCs in the outer region 
is consistent with the surface brightness profile of the 
stellar halo, (b) the elongation of the red GC system 
is consistent with that of the stellar halo, while the 
spatial distribution of the blue GC system is approxi- 
mately circular, (c) the mean color of the red GCs is 
similar to that of the stellar halo. These common prop- 
erties are also shown for other gE GC systems such as 
M87 and NGC 1399 (Bassino et al. 2006; Tamura et al. 
2006; Forte et al. 2007; Faifer et al. 2011; Strader et al. 
2011; Forte et al. 2012). Thus the red GCs are more 
correlated with the stellar halo than the blue GCs. This 
indicates that most red (metal-rich) GCs in gEs might 
have formed with the halo stars. 



5. SUMMARY 

We presented Washington CT\ photometry of the 
GCs in M86 covering a 16' x 16' field. We showed vari- 
ous photometric properties of the GC system. Primary 
results are summarized as follows. 

1. A significant population of GCs in M86 are shown 
in the CMD. 

2. The color distribution of the GCs in M86 is bi- 
modal, showing a blue peak at (C — T{) = 1.30 
and width a = 0.15, and a red peak at at 
(C — T\) — 1.72 and a = 0.15, with a division 
color (C - Ti) « 1.55. 

3. The spatial distribution of the blue GCs is roughly 
circular, while that of the red GCs is elongated 
similarly to that of the stellar halo. 

4. The radial number density profile of M86 GC sys- 
tem is more extended than the galaxy stellar light. 
The radial density profile of the red GCs is con- 
sistent with the surface brightness profile of M86 
stellar halo. The radial density profile of the blue 
GCs is more extended than that of the red GCs. 
The density profiles at 1' < R < 10' are approx- 
imately well fitted by a de Vaucouleurs law and 
a power law. As a result, the slope of the radial 
profile of the red GCs is steeper than that of the 
blue GCs. The radial density profiles in the cen- 
tral region are somewhat flat and approximately 
fitted by a King model. 

5. M86 GCs have the negative radial color gradi- 
ent because the number ratio of the blue GCs to 
the red GCs increases as galactocentric radius in- 
creases. The mean color of the red GCs is similar 
to that of the stellar halo. 

6. As seen in other Virgo gEs, the bright blue GCs in 
the outer region of M86 reveal a blue tilt that the 
mean colors of the blue GCs get redder as they 
get brighter. 
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